Pseudoscalar Goldstone bosons in the color-flavor locked phase at moderate densities 
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The properties of the pseudoscalar Goldstone bosons in the color-flavor locked phase at moderate 
densities are studied within a model of the Nambu-Jona-Lasinio type. The Goldstone bosons are 
constructed explicitly by solving the Bethe-Salpeter equation for quark-quark scattering in random 
phase approximation. Main focus of our investigations are (i) the weak decay constant in the chiral 
limit, (ii) the masses of the flavored (pseudo-) Goldstone bosons for non-zero but equal quark masses, 
(iii) their masses and effective chemical potentials for non-equal quark masses, and (iv) the onset 
of kaon condensation. We compare our results with the predictions of the low-energy effective field 
theory. The deviations from results obtained in the weak-coupling limit are discussed in detail. 



I. INTRODUCTION 

Much effort has recently been devoted to the study 
of strongly interacting matter at nonzero baryon den- 
sity. In particular the rich phase structure of color su- 
perconducting quark matter has attracted much inter- 
est. (For reviews on color superconductivity see, e.g., 
Refs. [H, H i, i, H i, 0, HO In nature, quark matter 
phases might be realized in compact stars [t| [13, EH- 
It is therfore natural to ask whether quark pairing has 
interesting phenomenological consequences for compact 
star physics. In this context the energetically lowest ly- 
ing degrees of freedom are relevant for many dynamic 
properties of quark matter. 

At low temperatures and very high densities the 
preferred state is most probably the color-flavor 
locked(CFL) phase where up, down, and strange quarks 
are paired in a particularly symmetric way [l2| . This can 
be shown from first principles within a weak-coupling ex- 
pansion [l3|, [TJ, [lj| . Although this expansion is not valid 
at "moderate" densities which could be reached in com- 
pact stars, recent Dyson- Schwinger studies indicate that 
the CFL phase might be the preferred phase all the way 
down to the hadronic phase [lq |. 

In the CFL phase, all quark flavors and colors partic- 
ipate in a condensate. As a consequence, all fermionic 
modes are gapped and do not appear in the low-energy 
excitation spectrum. The diquark condensates break the 
original U(l) baryon x SU(3) co i or x SU(3) L x SU(3) R sym- 
metry of three-flavor QCD (in the chiral limit) down to a 
residual Z-i x SU(3) co i or +Vi corresponding to a simul- 
taneous ("locked") rotation in color and flavor space. 
Due to the breaking of the color symmetry, all eight glu- 
ons receive a mass, while the breaking of baryon num- 
ber and chiral symmetry leads to the emergence of one 
scalar and eight pseudoscalar Goldstone bosons. In addi- 
tion, there is a ninth pseudoscalar Goldstone boson re- 
lated to the spontaneous breaking of Ua(1) which is a 
symmetry of QCD at very high density [l~7], [H[- In the 
presence of quark masses chiral symmetry is broken ex- 
plicitly and the pseudoscalar Goldstone bosons aquire a 
mass, while the scalar Goldstone boson remains mass- 



less. Since, with all quarks being gapped, the Goldstone 
bosons are the lowest lying excitations, they play an im- 
portant role for the thermodynamic and transport prop- 
erties of strongly interacting matter, relevant for compact 
star phenomenology (cf., e.g., Refs. fl9l. I20I [2lT|). 

The symmetry breaking pattern is the basis for 
the construction of the low-energy effective theory 
(LEET) [H, HI HI HE H, S3, l2|, which describes 
the Goldstone boson dynamics and is valid for ener- 
gies much smaller than the superconducting gap. At 
very high densities, the interaction is weak and the con- 
stants for the LEET can be calculated from QCD using 
High Density Effective Theory (HDET) [H, El H EI . 
For instance, in the weak-coupling limit, pseudoscalar 
meson masses and deca y c onstants have been investi- 
gated [H, H3, H3, [H, [H, [111 . It was also shown that the 
stress imposed by the strange quark mass on the CFL 
Cooper pairs acts as an effective strangeness chemical 
potential, which may eventually lead to kaon condensa- 
tion [H,[2S[13- 

At intermediate densities, relevant for compact 
star phenomenology, the interaction becomes non- 
perturbative and it is difficult to study the Goldstone bo- 
son dynamics from first principles. The leading-order pre- 
dictions, however, are often universal, in the sense that 
they do not depend on the interaction, but should hold in 
any model exhibiting the same symmetry pattern. One 
such model is the Nambu-Jona-Lasinio (NJL) model [36| . 
often used in the intermediate density regime to study at 
least qualitatively the main features. (For reviews, see, 
e.g., Refs. d, [I?], [H, HU.) This model has already been 
applied to study kaon condensation in the CFL phase 
at non-zero strange quark masses [4(| [U [42[ . However, 
this was done by focusing on the ground state proper- 
ties, i.e., without explicit construction of the Goldstone 
bosons. In Refs. [H, [44|, on the other hand, meson and 
diquark properties in the CFL phase have been studied 
explicitly, but this investigation was restricted to the chi- 
ral limit. (Mesons and diquarks in the 2SC phase have 
been discussed in Refs. [45l.l46|.) 

In the present paper we discuss a detailed analysis of 
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properties of pseudoscalar mesons 1 in the CFL phase in 
an NJL-type model including the cases of equal and un- 
equal nonzero quark masses. Emphasis is put on a com- 
parison with the weak-coupling results. 

The paper is organized as follows. In Sec. [IT] we intro- 
duce our model and discuss how to construct the mesonic 
excitations. Sec. IIIII is devoted to some general results 
which can be obtained in the limit of equal quark masses 
based on chiral Ward-Takahashi identities. In Sec. lIVI nu- 
merical results will be presented. Within that section, we 
investigate the pion decay constant in the chiral limit as 
well as meson masses for the cases of equal and unequal 
quark masses. In this context we also discuss the onset 
of kaon condensation in the CFL phase. Our results are 
summarized in Sec. [V] 

II. FORMALISM 

A. Model Lagrangian 

We consider an NJL-type Lagrangian 

£ = q{ift- m)q + C qq , (1) 

where q is a quark field with three flavor and three color 
degrees of freedom, rh = diagj(m„, mj, m s ) is the mass 
matrix, and 

£qq=H 2_j (qij5TA\A>Cq T )(q T Cij 5 TAXA>q) 

A,A'=2,5,7 

+ (qTA\A'Cq T )(q T CT A \A>q)] (2) 

describes an SU(3) co ior X U(3)l X U(3)r symmetric four- 
point interaction with a dimensionful coupling constant 
H. C = ij 2 j° is the matrix of charge conjugation, and r 
and A, denote Gell-Mann matrices acting in flavor space 
and color space, respectively. In this article, we follow 
the convention that the indices A and A' are used for 
the antisymmetric Gell-Mann matrices only, i.e., A, A' £ 
{2,5,7}, whereas arbitrary Gell-Mann matrices will be 
denoted by small Latin letters, e.g., r a , a = 1, . . . , 8. 

The first term in Eq. ^ corresponds to a scalar quark- 
quark interaction in the color and flavor antitriplet chan- 
nel, just as needed for giving rise to the diquark conden- 
sates in the CFL phase (see Eq. (fT2"|) below) . The second 
term is the corresponding pseudoscalar interaction and 
is required by chiral symmetry. This will be the essential 
term for the pseudoscalar Goldstone excitations we want 
to study. 

For simplicity, we restrict ourselves to quark-quark 
interactions. The effect of quark-antiquark interactions, 
which give rise to normal self-energies and thereby to 
dynamical quark masses, will be investigated in a future 
publication. 



1 In this article, we often use the words "mesons" and "diquarks' 
synonymously, see Sec. Ill E51 for more details. 



B. Operators in Nambu-Gorkov space 

Introducing Nambu-Gorkov bispinors, 

* - (<*) ■ < 3 > 

Eq. ^ can be rewritten as 

c qq = ah { *r^**r^,* 

A,A'=2,5,7 

+ *r£ 1 ,**r£ 4 ,*} 1 (4) 

with 18 scalar operators 

r sT _ (0 «75TaAaA r si _ ( 0\ (v s 
1 aa> - [q o )> l aa> ~ \i l5TA X A , o) W 

and 18 pseudoscalar operators 

r pt _ (° t a X a \ r pl . / Q 0\ , , 
1aa ' ~ [p J' 1 aa- ~ \r A \ A ,ty- W 

From these expressions we obtain the scattering kernel 

K = TiKi^j, (7) 

where Tj are the 36 operators defined above, 

f, - 7°rl 7 , (8) 

and 

Ka = AH5 l3 . (9) 

Repeated operator indices are summed over, unless 
stated otherwise. 

Vertices describing the coupling of an external source 
to a bare quark are generalized to Nambu-Gorkov space 
in the following way: 

f ^»« s 6- c ?4 < io > 

This guarantees that quark-antiquark bilinears remain 
unchanged, (f) atg ^ = <?r <?• 

C. CFL ground state 

Before we construct the mesonic excitations, we have 
to determine the ground state of the system. For equal 
quark masses the CFL phase can be characterized by the 
equality of three scalar diquark condensates in the color 
and flavor antitriplet channel, 

S22 = S55 = S77, (11) 

where 

SAA' = (q T Cj 5 t a \a> q) ■ (12) 



In general, these condensates are accompanied by in- 
duced color-flavor sextet condensates. These are, how- 
ever, small and can be neglected. If the SU (3)-flavor 
symmetry is explicitly broken by unequal quark masses, 
Eq. (fTTj) does no longer hold exactly, but the three con- 
densates may differ from each other. 

To obtain the ground state we must minimize the ther- 
modynamic potential (per volume V), 



T 
V 



In Z(T, {/*}), 



(13) 



where Z(T, {/li}) is the grand partition function at tem- 
perature T and a given set of chemical potentials {fii}- 
For P equilibrated matter, these can be expressed in 
terms of the quark number chemical potential /j,, the elec- 
tric charge chemical potential fiQ, and two color chemical 
potentials ^3 and n$ (47| . 

For the CFL phase in mean-field approximation, O is 
given by 



+ 177 E l A ^l 2 ' ( 14 ) 



A=2,5,7 



where the gap parameters A^ are related to the diquark 
condensates, 



A A = -2Hs AA . 
The inverse dressed propagator reads 

i> + At - Z) ^ai^aX, 



(15) 



s-\v) = 



4=2,5,7 

A^taXa f - /t7° - to 

VA=2,5,7 



(16) 

Here /x denotes the diagonal matrix in color and flavor 
space which is given by the set of chemical potentials 

{«}• 

In order to determine the ground state, f2 must be min- 
imized with respect to the gap parameters A a, leading 
to three gap equations 



on 



dA 



- = 0, A = 2,5,7. 



(17) 



Furthermore, we require the solutions to be color and 
electrically neutral in the presence of leptons. In general, 
this leads to three additional equations, 



= 0, i = Q, 3,8, 



(18) 



where VL to t = £2 + ^leptons is the sum of the quark part, 
Eq. (fT4)) . and the contribution of the leptons. Thus, al- 
together we have a set of six coupled equations for A^ 
and ^ which must be solved simultaneously. This has 
been done many times before, and we can refer to the 
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FIG. 1: Dyson equation for the dressed Nambu-Gorkov quark 
propagator (thick line). The thin line indicates to the bare 
propagator. 



literature for technical details, e.g., Refs. [471. l48l. I49u50l] . 
In the present article, we restrict ourselves to the (fully 
gapped) CFL phase at zero temperature. In this case, 
color neutral quark matter is electrically neutral without 
leptons (5lj and /iq = 0. Moreover, we consider isospin 
symmetry, m u = m c i, so that 113 vanishes as well. 

The gap equations (|17[) can be derived from the Dyson 
equation for the dressed quark propagator, too, diagram- 
matically shown in Fig. [TJ This is well known, but some 
details are useful in our later discussion. Therefore, we 
present this derivation in Appendix |A"1 



D. Axial transformations 

In the chiral limit, m u = to^ = m s = 0, the La- 
grangian, Eq. |T]), is invariant under SU (3) co i or x U(3)l X 
U(3)r transformations. In the CFL phase, this symmetry 
is spontaneously broken to the diagonal vector subgroup, 
SU(3) co i or +y . This is reflected by nonzero values for the 
condensates, Eq. (jlip . Axial transformations, 



(I 



exp(i6> a 7 5 i a )g, 



(19) 



then connect a continuous set of degenerate ground states 
for the CFL phase. Here t a = a — 0, ... ,8, where 
T\ , . . . , Tg denote the eight G cll-M ann matrices in flavor 
space, as before, and 7o — \/2/3 1 /. These transforma- 
tions can be parameterized by nine pseudoscalar Gold- 
stone bosons. In order to identify the operators corre- 
sponding to Goldstone boson excitations, we inspect the 
effect of an axial transformation with specified quantum 
numbers on the condensates, Eq. (|lip . In Nambu-Gorkov 
formalism this can be described as follows 

(§Fi*) -> (*'ri*') = (*r<*), (20) 

for Ti = T S 2 A and T s j A . Conside ring infinitesimal trans- 
formations, we find 

^ = r< + ie a sr i>a , (21) 

with 

ST^ = {(76t )jVG, Ti}, (22) 
where {A, B} = AB + BA denotes the anticommutator. 
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TABLE I: Pseudoscalar Goldstone modes, corresponding fla- 
vor operators tj, and Nambu-Gorkov operators Tj contribut- 
ing to the quark-meson vertex. The Tj are obtained by eval- 
uating Eq. (f22)l for t a = tj and all possible T; = T"2 A or T AA . 

x = x+xx 

FIG. 2: RPA equation for the T-matrix in Nambu-Gorkov 
space. 



The Goldstone bosons are coupled to the quarks by 
vertices related to <5I\ a (see Ref. [32[ for a detailed dis- 
cussion). Thus, by evaluating Eq. (f2"2")l for I\ = T"J A and 

we can determine the operators which contribute 
to the vertex of a given Goldstone boson. It turns out 
that the assignment is most simple if we work in the 
"particle basis" instead of using hermitian flavor opera- 
tors. In this case all flavored mesons are coupled to only 
two Nambu-Gorkov operators, while the hidden-flavor 
mesons are coupled to four (n°) or six (r}$,r}o) operators. 
The results are summarized in Table [I] As we will see be- 
low, this assignment remains basically unchanged when 
we include finite quark masses. 



E. Mesonic excitations 

Iterating the scattering kernel as illustrated in Fig. [2] 
leads to the RPA equation for the T-matrix in Nambu- 
Gorkov space, 

f = K + KJ f. (23) 

Since the interaction is separable, this operator equation 
can be converted into a matrix equation by using Eq. ([7]) 
for the scattering kernel K, an analogous expression for 
f, 

f = \ ,T,A r (24) 



and defining 

ViJTj = Jij. (25) 

One finds 

T = K + KJ T, (26) 

with the solution 

T(q) = (l-KJiq))- 1 K = (± -J(q)^ , (27) 

where the last equality follows from Eq. (|9|). The polar- 
ization matrix, corresponding to the loop in Fig. [2l is 
given by 

JiM = * / \^ [ f iS{k + q)T 3 S(k)] . (28) 

Here we have introduced a "vacuum-like" notation for 
brevity. In medium we should replace 



and 




with bosonic Matsubara frequencies ioj m and fermionic 
Matsubara frequencies iio n . In the end, the result should 
be analytically continued to real external energies. We 
will use this notation throughout this paper. For our nu- 
merical results, we will introduce a 3-momentum cutoff 
A to regularize the divergent integrals. 

The matrices T and J are 36 x 36 matrices in opera- 
tor space, corresponding to the 36 operators defined in 
Eqs. ((5|) and ((6]). It turns out, however, that J and, thus, 
T are block diagonal, i.e., only certain combinations of 
operators occur. More precisely, scalar operators do not 
mix with pseudoscalar ones, and each of the resulting 
18 x 18 blocks can be decomposed further into six 2x2 
blocks and one 6x6 block. These blocks carry different 
quantum numbers and reflect the assignment to diffe- 
rent meson modes, as given in Table |T] 2 In particular, 
the 2x2 blocks correspond to the flavored mesons. The 
6x6 blocks, on the other hand, describe the hidden-flavor 
mesons, i.e., 7r°, r/% and 770 in the pseudoscalar case. As in 
vacuum, these states are mixed for unequal quark masses. 

This assignment can be tested by coupling the T- 
matrix to an external meson source, as illustrated in 
Fig. [3] To that end, we evaluate the loop integral 



2 In Table U we have listed the pseudoscalar mesons only, but the 
scalar sector is completely analogous. 
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FIG. 3: Coupling the T-matrix (double line) to an external 
meson source (wavy line). 



hj{q) = i 



d 4 k 1 
(2^*2 



:S(k + q ) (rr) NG S(k) 



(31) 

where Tj xt denotes the vertex of the external source to 
the quark. It is given by Tj xt = tj for scalar sources 
and Tj Xt = ij 5 tj for pseudoscalar sources. Here tj is one 
of the generators in flavor space listed in Table [I] and 
(r^ x *)jvG is the generalization of this vertex to Nambu- 



Gorkov space, as defined in Eq. (flO)) . 

The blocks can be diagonalized by unitary transforma- 
tions, 



J' = WJW\ 



(32) 



with W being an (in general 4-momentum dependent) 
unitary matrix and 



Defining new operators 

J' can be rewritten as 

d 4 k 1 



(33) 



(34) 



JIM 



(2tt) 4 2 



Tr [T'iSik + g)rj5(fc)] . (35) 



Since the scattering kernel remains diagonal in this new 
basis, 



k = 4iir,r 4 = 4ijr^n, 

the T-matrix becomes diagonal as well 
with 



i 

iff 



(36) 



(37) 



(38) 



In the vicinity of a pole, we can parameterize these modes 
like a free boson with mass m; in the presence of a chemi- 
cal potential //$ corresponding to this particular boson, 



"5i 



(go + ^) 2 -cfg 2 



(39) 



Here <?i can be interpreted as a coupling constant of the 
boson to an external quark, and Cj denotes the in- medium 



group velocity. In this article we restrict ourselves to q = 
in order to keep artifacts of the 3-momentum cutoff as 
small as possible. 

The various modes T^ 1 ' describe bosonic excitations 
of the CFL ground state. Because of the formal quark- 
antiquark structure of the polarization loops, we will call 
these excitations "mesons". However, it should be kept 
in mind that the propagators and vertices entering the 
polarization loops live in Nambu-Gorkov space and there- 
fore in principle describe quark-antiquark as well as di- 
quark and antidiquark (or di-hole) excitations. In vacuum 
or in a normal-conducting medium, these are indepen- 
dent modes, protected by the conserved baryon number. 
In the CFL phase, however, baryon number is broken, 
and quark-antiquark, diquark, and antidiquark states can 
mix. As our model Lagrangian does not contain quark- 
antiquark interactions, our "mesons" are in fact superpo- 
sitions of diquarks and antidiquarks or, more important, 
di- holes. 

In vacuum we have nine scalar and nine pseudoscalar 
diquarks and nine scalar and nine pseudoscalar antidi- 
quarks. Since the total number of states does not change 
when the states are mixed, there must be 36 meson states 
in the CFL phase, 18 scalars and 18 pseudoscalars. (If we 
included quark-antiquark interactions in our model, we 
would obtain 27 scalars and 27 pseudoscalars.) As we will 
see below, nine pseudoscalars are massless in the chiral 
limit, while the others stay heavy. This has been found 
in Ref. [44|, too, within a similar model. 



F. Pseudoscalar meson decay constants 

As in vacuum, the pseudoscalar mesons can decay 
weakly. The decay amplitudes are related to the loop in- 
tegrals 



(a) 



d A k 1 



(2tt) 4 2 



Tr^ S(k + q) (ri5tj)N G S(k)) , 

(40) 

describing the coupling of the meson i (i.e., the one which 
corresponds to the operator T^) to an external axial cur- 
rent A^j. 3 In the following, we are mostly interested in 
flavored mesons, which are the main focus of our studies. 
Each flavored meson i couples to only one tj, as listed 
in Table HI 4 For simplicity, we denote the flavor operator 
which fits to the meson i by U. 



3 Strictly speaking, Eq. (1406 describes the time-reversed process, 
i.e., the production of a meson by an incoming axial current. This 
choice was made for later convenience in Sec. IIII Bl There we 
apply chiral Ward identities, which are formulated for incoming 
axial currents, see Eq. J52J. 

4 Note that the opposite is not true: Since we have 18 pseu- 
doscalars, there are in general two meson states i which couple 
to a given flavor operator tj . 



The meson decay constants are related to the on-shell 
values of these amplitudes, 



on — shell 



(41) 



with no summation over the index i on the r.h.s. The cou- 
pling constant gi has been defined implicitly in Eq. (|39|) . 

In general, there are different values for the time-like 
(/j = 0) and the space-like (// = 1,2,3) decay con- 
stants, which differ by the group velocity of the Goldstone 
modes. However, since we only consider mesons with van- 
ishing 3-momenta in this article, we are restricted to the 
time-like decay constants. 



III. EQUAL QUARK MASSES 









1 ; E 




FIG. 4: Vertex function for an external axial current (wavy 
line) coupled to a quark. The shaded circles (black dots) in- 
dicate dressed (bare) vertices. The open circles correspond to 
quark-quark vertices and the double line to the T-matrix, as 
in Fig. [2] 



It is rather instructive to investigate the simplified case 
of equal quark masses, m u — rrid = m s = m. In this case, 
we have only one gap parameter A 2 = A 5 = A7 = A and 
one common chemical potential /i in electrically and color 
neutral CFL matter. Moreover, the set of the 18 pseu- 
doscalar meson states consists of two 5/7(3) octets and 
two singlets, with all mesons in the same multiplet being 
degenerate. 5 Finally, there is no stress caused by quark 
mass or chemical potential differences, which could act 
as an effective meson chemical potential as in Eq. ([39| . 
Hence, for vanishing 3-momenta, the pole approximation 
for the T-matrix, Eq. (f3"9")l . becomes 



(42) 



T {t) (qo) 



1 2 o 



(43) 



This yields for the decay constants 

qa=mi,q=Q 

We will often refer to the states of the lowest pseudoscalar 
octet as "pions" and denote their masses, couplings, and 
decay constants by m^, g n , and f^, respectively. 

A. Dressed vertex functions 

Attaching an external axial current to both sides of 
the Dyson equation for the dressed quark propagator (see 
upper line of Fig. [J), we obtain a sclfconsistency equation 
for the dressed axial vertex, 



+ AiHTi 



(7^75 ij)wc 

Al I lTr[f i S(k + q)T^(k;q)S(k) 



5 The degeneracy in the octets is due to the residual 5(7(3) color+v 
symmetry of the CFL phase. Thus, even though we start from 
a U(3) invariant Lagrangian in our model, there is no unbroken 
symmetry which relates the singlet states to the octets. 



The first term on the r.h.s. corresponds to the bare ver- 
tex, while the second term contains the dressed vertex 
again. Thus, iterating the equation, the dressed vertex 
can be written as a Born series of quark-antiquark loops 
(second line of Fig. Q| , with the bare vertex attached to 
the last loop. Employing the quark-antiquark T-matrix, 
this could be rewritten (last line of Fig. 0|) as 

rfcfag) = (i^Wng + T«(g)r^/(«), (45) 

with as defined in Eq. (|4"0")) . We have seen earlier 
that a given flavor operator tj couples to only two meson 
modes i in flavored channels. We may call these modes 
"•7Tj" and ££ 7r'-". Then, for q = and go approaching m n , 
we obtain 

r?,-(p; qo — * m Tr) — i ^j— - r' + non-singular terms. 
q£ - mi 1 

(46) 

Similarly, we have an equation for the dressed pseu- 
doscalar vertex 

r 5j(p;<7) = (ij5t 3 )NG 

d 4 k 1 



4iH r, 



r) i4 -Tr[T i S(k + q)T 5j (k;q)S(k)] , 



(47) 



which could be rewritten as 

T 5j (p;q) = (i l5 tj) NG + T (l \q)r' l I' l] (q), (48) 

with 

w f , . r d A k 1. 



(2tt) 4 2 



Tr [T'iSik + q) (i^t^NoSik)] 



(44) 

Thus, approaching the pole, we find 



(49) 



2 j(7r) 



^5j(p;qo -* rn n ) = — I 1 

qo-ml 



with 



r' . + non-singular terms, 
(50) 



(51) 
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B. Chiral Ward-Takahashi identity 



for the flavored mesons, and 



It can be shown on general grounds that the exact 
vertex functions and inverse propagators must satisfy the 
following axial Ward-Takahashi identity: 



= S^ 1 (p + q) (75 *j)tfG + (l5tj)NGS~ 



L (P)- (52) 



In Appendix [Bl we show by using the gap equation that 
the vertex functions defined above and the inverse prop- 
agator (fl6|) are consistent with this relation. 

We now evaluate this equation for q = and qo — > 
m v . In the chiral limit, m = 0, only the axial vertex 
contributes, and we find from Eq. ([46)) that the l.h.s. is 
given by 



% 



5j 
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g^f-n^n + non-singular terms. 



In general, this has a singularity at go = m -n ■ On the other 
hand, one can easily see that the r.h.s. of Eq. ([52)) remains 
finite. We thus conclude that the singularity on the l.h.s. 
must be suppressed, i.e., either m t = or g v f n = 0. 
In fact, from the symmetry breaking pattern we expect 
nine pseudoscalar Goldstone bosons. Therefore, both sce- 
narios should be realized, i.e., nine pseudoscalar mesons 
(one octet and one singlet) are massless and the others 
are massive. 

For the massless solution, we can now evaluate Eq. 
directly at q = 0. This yields 



(54) 



Inserting Eq. (|16|) for the inverse propagator with m = 
and taking the freedom to choose the gap parameter A 
to be real, we find the solution 
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A=2,5,7 



(55) 
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)}• (56) 



Eq. (|55[) may be viewed as a generalization of the well- 
known Goldberger-Treiman relation in vacuum (32j. 
Explicit evaluation of Eq. ([56)) yields 
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for the mesons with hidden flavor, 
modes are related to each other as 

p'_ (o) = f / (0) 



(58) 

Thus, antiparticle 
(59) 



(with 7T°, 778 and 770 being their own antiparticles) as it 
should be. 

The results above are consistent with the operators 
identified in Ref. [32| following a similar reasoning. Ex- 
cept for an arbitrary phase, Eqs. ([57)) and ([58")) also agree 
with the result of diagonalizing the meson polarization 
function Jij for vanishing quark masses (see Sect. Ill El) . 
We should keep in mind, however, that we have evalu- 
ated the meson vertices at the pole only. In general, the 
diagonalization of the T-matrix leads to 4-momentum de- 
pendent vertex functions with 4-momentum dependent 
weights for the contributing operators. For instance, the 
7r + vertex can be written as 



r' 



(q) = -t(Bin^(g)rg| - cos <p(q)I%i), (60) 



with some function <p(q) and an arbitrary overall phase 
which we have chosen to be purely imaginary to comply 
with Eq. (|57p. For equal quark masses, <p(0) = \ and we 
recover 1^(0) = T'J® , Eq. ([57]). 

Next, we evaluate Eq. ([52)) for non-vanishing (but still 
equal) quark masses m. In this case the pseudoscalar ver- 
tex contributes a pole at q = and qo — > m n , which is 
present for any choice of (see Eq. ([50)) ). Since the 
r.h.s. of Eq. ([52")) remains non-singular, this pole must 
be cancelled by the pole in the axial part (unless g 2 1^' 
vanishes) . This means that the formerly massless mesons 
receive a mass which is determined by the requirement 
that the residues of the axial and the pseudoscalar pole 
cancel each other. One finds 



ilfl = 2m 9 ,/jW. 



(61) 



This is an exact relation, valid for arbitrary values of m. 
We can now perform a chiral expansion of this equation 
to leading order. This amounts to replace g^f^ by A and 
to evaluate 
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-Tr[fUg)S(fc- 
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to leading nontrivial order in the quark mass. It turns 
out that jW vanishes in the chiral limit and the lead- 
ing order is linear in m. 6 Note that, in principle, /W 
depends explicitly (via the quark propagators) and im- 
plicitly (via m w ) on m. One can show, however, that the 
implicit contributions vanish in leading order. It is there- 
fore consistent to evaluate the integral at go = 0, and we 
obtain 



with 



A = — i 
4 



■ J TV 



! Am + higher orders, 



(63) 
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(64) 



and r^°) = T^(O) as given in Eq. (|57|) . This expression 
can be evaluated exactly. The result is given in Eq. (|E1[) 
in the appendix. Expanding that formula in A, we find 



.4 



A 



^2 (-hit/ 2 -2 + ^ln2 + ln(a; 2 -l)) + 



where we have introduced the abbreviations 



A 

x = —. 



A 

y = — , 



(65) 



(66) 



with A being the 3-momentum cutoff. This should be 
compared with the weak-coupling result [23[ , 



A, 



3A 2 
4^" 



(67) 



An important difference is the logarithmic term in 
Eq. (f6"5)) . which does not exist in Eq. (f6"7) . In fact, in 
the beginning, there was a controversy about the correct 
weak-coupling limit, and similar logarithmic terms have 
also been discussed in the literature [H, H3, [H, [H, [H], 
[52| . To be precise, it was found that the leading contri- 
bution to A in weak coupling is given by [30l ] 
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AA 
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(68) 



where A is the antiparticle gap. In the NJL model, the 
gap function is energy independent and, hence, A = 
A. Therefore, our result, Eq. (|65|) . is consistent with 
Eq. (|68p. Moreover, if we introduce particle and antipar- 
ticle gaps by hand as independent constants in the quark 
propagator, we recover Eq. ([68]) in leading order. 



If we had taken into account dynamical quark masses it would 
be linear in the constituent quark mass M. The leading term on 
the r.h.s. of Eq. J63H would then read 8AmM. 
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FIG. 5: Pion decay constant f n in the chiral limit as a func- 
tion of the gap parameter A: numerical results (points) in 
comparison with the semi-analytical formula Eq. (|70|l . ne- 
glecting (thin lines) and including (thick lines) the momentum 
dependence of the vertex function. Asterisks and solid lines: 
A = 600 MeV; crosses and dashed lines: A = 700 MeV; dotted 
line: weak-coupling limit, Eq. (1691 1, 



On the other hand, as argued in the erratum of Ref. [2j| 
and confirmed in Ref. [53| , in a gauge invariant treatment 
of QCD at weak coupling, the antiparticle gap contribu- 
tions are cancelled by other terms, and the logarithmic 
terms drop out. This finally leads to Eq. (|6"T|) . Of course, 
gauge invariance is not an issue in the NJL model. Nev- 
ertheless, the logarithmic behavior at very weak coupling 
should be viewed as a model artifact. 



IV. NUMERICAL RESULTS 



In the following we present our results for T 
a fixed quark chemical potential \i = 500 MeV. 



and 



A. Pion decay constant in the chiral limit 

We begin with a discussion of the pion decay con- 
stant in the chiral limit, i.e., for vanishing quark masses. 
In Fig. [5] our results are displayed as functions of the 
gap parameter A for two different values of the cutoff 
A. The points correspond to the numerical results for 
A = 600 MeV (asterisks) and A = 700 MeV (crosses). To 
be precise, these calculations have been performed with 
m u — rrid = m s =0.1 MeV for technical reasons. 

The weak-coupling limit, A — * 0, of f n has been de- 
rived in Ref. [23| from an effective theory involving only 
fermionic modes in the vicinity of the Fermi surface. The 
result, 
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(69) 
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is universal and should hold in any model exhibiting the 
same symmetry pattern. 7 Indeed, our results converge 
to this limit for A — > 0. In the general case we find, 
however, deviations from Eq. (f6"9")) . e.g., about 10% for 
A = 600 MeV and A w 80 MeV. Moreover, these devia- 
tions depend rather strongly on the cutoff. 

In order to understand this behavior and to confirm 
the correct weak-coupling limit, we employ the semi- 
analytical formula derived in App. |Dl 



fl = fl + Sf 2 . 



(70) 



f 2 thereby describes the contribution to /„. arising for a 
constant pion vertex function, i.e., neglecting the energy 
dependence of the mixing angle ip in Eq. (|60[) . This part is 
given in a closed analytical form in Eq. (|D5[) . Sf 2 , on the 
other hand, incorporates the effect of the energy depen- 
dence of the vertex function and is given in Eq. (|D1 1|) . It 
is proportional to the derivative dip/dqo, which is evalu- 
ated numerically. 

Expanding / 2 and the analytical factor of Sf 2 for small 
values of A yields 



f 2 = {(21 -8 In 2) - 9y 2 lny 2 - c 2 y 2 + . .. } (71) 



and 



Sfl = -^{-36yW-4 
where 
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dp 
dq 
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> (72) 



c, = 181n2-9m(a; 2 -l) 
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A5x 2 - 27 



(x 2 - l) 2 ' 

d 2 = 102 - 561n2- 361n(x 2 - 1) + ^ , (73) 

ar — 1 

and x and y are defined in Eq. (f6l>)) . 

We see that the weak-coupling limit, Eq. (|69|) . is cor- 
rectly reproduced by the leading term in / 2 , while Sf 2 
does not contribute to this order, provided dp/dqo does 
not diverge as strongly as {y 2 h\y 2 )~ 1 for A — » 0. We 
therefore expect that f 2 gives the main contribution to 

As already mentioned, the weak-coupling limit is uni- 
versal and therefore must be cutoff independent. In the 
detailed calculations, this results from the fact that in 
the limit A — > the integrand of the 3-momentum in- 
tegral in Eq. (ID3[) becomes proportional to <5-functions 
at the Fermi surface. While the y 2 lny 2 -term in Eq. (|7Tj) 
is cutoff independent as well, the quadratic term is not 
universal and dependends on A via the variable x. This 



7 Strictly speaking, the chemical potential {i should be replaced by 
the Fermi momentum pp because in general the Fermi velocity 
can differ from the speed of light. 
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FIG. 6: Derivative of the mixing angle tp as a function of the 
gap parameter A. Solid: A = 600 MeV, dashed: A = 700 MeV. 



is also the case for the quadratic term in Eq. (|T2[) . Any- 
way, the situation for Sf 2 is more complicated because 
the derivative dip/dqo depends on the cutoff as well. 

The results of the semi-analytical formula for /„. are 
indicated by the thick lines in Fig. \5\ The solid line 
corresponds to A = 600 MeV and the dashed line to 
A = 700 MeV. We have employed the exact formulas, 
Eqs. (|E>5|) and (jDlip . for f 2 and Sf 2 , respectively, with 
dip/dqo being computed numerically. Obviously, the re- 
sults for /„• are in perfect agreement with the numerical 
computations. 

In order to analyze the influence of the momentum 
dependence of the vertex function, we display the func- 
tion /tt, too (thin lines). Since f w contains the leading 
term, Eq. (|69p . it is not surprising that it is the domi- 
nant contribution to However, for a correct descrip- 
tion of the deviations from the weak-coupling limit, 5f 2 
can be quite important: Whereas for A = 600 MeV, we 
find that f n is well reproduced by / T , this is not the 
case for A = 700 MeV. This indicates that the correc- 
tion due to the contributions arising from the momentum 
dependence of the vertex function are rather small for 
A = 600 MeV and considerably larger for A = 700 MeV. 

This result becomes clear if we look at the deriva- 
tive dip/dqo at q = 0, which is shown in Fig. [5] For 
A = 700 MeV (dashed line) the derivative is nowhere 
small in the shown region. Therefore, Sf 2 can never be 
neglected, in agreement with our findings in Fig. [SJ On 
the other hand, for A = 600 MeV (solid line) we find that 
the derivative is rather small for A > 40 MeV, explaining 
why Sf 2 is negligible in this regime. For smaller values 
of A the correction becomes larger, in agreement with 
the deviations we found in Fig. In fact, for A — > 0, 
the derivative even seems to diverge. If it was stronger 
than (y 2 In J/ 2 ) -1 , the divergence could affect the weak- 
coupling limit. However, at least numerically, we find 
that dp/dqo grows much slower than (y 2 lny 2 ) -1 , and we 
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FIG. 7: fir (solid line) and A (dashed line) divided by their 
respective chiral limit values and A' ' as functions of the 
common quark mass m. The calculations have been performed 
for HA 2 = 1.4, corresponding to A (0) = 79.1 MeV. 
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FIG. 8: Masses of the flavored (pseudo-) Goldstone bosons as 
functions of a common quark mass m for three different cou- 
plings. Solid: HA 2 = 1.4 (corresponding to A = 79.1 MeV), 
dashed: HA 2 = 1.0 (A = 43.2 MeV), dotted: HA 2 = 0.6 
(A = 12.5 MeV). 



therefore conclude that the weak-coupling limit is safe. 



B. Equal quark masses 

Next, we study the effect of explicit chiral symmetry 
breaking through non- vanishing, but equal quark masses 
m u = rrid = m s = m. For the cutoff, we choose A = 
600 MeV and keep this value fixed in the remainder of 
this article. 

Before turning to our main focus, i.e., the masses of the 
Goldstone bosons, we briefly investigate the dependence 
of the pion decay constant and of the gap parameter A on 
m. The results are displayed in Fig. \7\ We find that both 
quantities depend only very weakly on the quark mass. In 
the plotted range they vary less than 0.2 %. This is much 
weaker than the m dependence of / ff in vacuum or of the 
constituent quark mass in comparable models. In the fol- 
lowing discussion we will therefore neglect the distinction 
between the m dependent gap parameter A and its chiral 
limit A' ' and often use A in order to characterize the 
coupling strength. 

From Eq. (f6"5)) . we expect that, to leading order, the 
masses mjf of the Goldstone bosons in the octet depend 
linearly on m, 



m M 




(74) 



This is confirmed by our numerical calculations. In Fig. [5] 
the values of mu are displayed as functions of the quark 
mass for three different couplings HA 2 = 0.6, 1.0, and 
1.4, corresponding to CFL gaps of 12.5 MeV, 43.2 MeV, 
and 79.1 MeV, respectively. As one can see, the results 
show an almost perfect linear behavior. 



The slopes of the straight lines are decreasing with 
decreasing coupling strength H, i.e., with decreasing A. 
This is also expected from the QCD weak-coupling limit, 
Eq. (jnZl) and Eq. |g| inserted in Eq. (J7U). However, as 
discussed in Sect. lIIIBl we expect that the A-dependence 
of the slopes in the NJL model and in weak-coupling 
QCD is rather different. To study this aspect in detail, we 
determine the slopes a of the functions rriM^jn) = am for 
different values of A and use our (A dependent) results 
for to obtain A = ^a 2 f 2 . 

The result, divided by the weak-coupling limit A wc is 
displayed in Fig. [5] as a function of A. The ratios ob- 
tained from the fit to the numerically determined meson 
masses are indicated by the crosses. The solid line corre- 
sponds to Eq. (|E1|) , i.e., to the exact analytical solution of 
Eq. (|64[) . the dashed line indicates the approximate for- 
mula Eq. (|65[) . The former is again in perfect agreement 
with the numerical results. 

We see that for small values of A, A is larger than A wc . 
This is due to the logarithmic term discussed in Sec. lIII Bl 
(see Eq. (1651) and the subsequent discussion). For large 
couplings, on the other hand, the NJL-model value of A 
is considerably smaller than A WC1 leading to even smaller 
Goldstone-boson masses than predicted in weak coupling. 
This has also been found in Ref. [4(| within a rather 
different approach. 



C. Unequal quark masses 

Finally, we study the effect of unequal quark masses. 
In the upper panel of Fig. [10] the pole positions of the 
flavored (pseudo-) Goldstone modes are displayed as 
functions of the strange quark mass, keeping m u and 
rrid fixed. Our numerical results are indicated by the 
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FIG. 9: Ratio of A and A wc as a function of the gap parameter 
A. Points: numerical results; solid line: exact analytical result 
with A from Eq. (|E1[) ; dashed line: approximate result with 
A from Eq. (|65) . 



points. For practical reasons, we have chosen a relatively 
strong diquark coupling HA 2 = 1.4 (corresponding to 
A = 79.1 MeV) and a relatively large value of 30 MeV for 
m u and mj, in order to have not too small meson masses. 
As one can see, the poles which are degenerate for equal 
masses split into three branches, corresponding to diffe- 
rent strangeness, i.e., pions (5 = 0), kaons (S = +1), and 
antikaons (S = —1). On the other hand, since we have 
chosen m u — m^, the different isospin states of these 
modes, i.e., ir~ and ir + , K + and K°, and K~ and K° 
remain degenerate. 

Our results can be analyzed in terms of the parameter- 
ization given in Eq. (|3"9"|) . Each mode TW has two poles, 
which for q = are located at 



qo = ±m,i 



(75) 



We can thus extract the meson masses and chemical po- 
tentials as 



mi = -« - ^ ), 



1 



(ut + u>r). (76) 



The resulting functions are displayed in the two lower 
panels of Fig. 1101 The masses are moderately rising with 
m s and exhibit an "inverse ordering" (mx < m^), as 
predicted first in Ref. [23|. Also note that kaons and an- 
tikaons have the same masses. For the meson chemical 
potentials, on the other hand, we find fj,K — as it 

should be, and fj,^ = 0. 

These results can be compared with those derived in 
Refs. [23] and [H, H3 in an effective field theory (EFT) 



approach. For m u = m<j 



they read 
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FIG. 10: Properties of the flavored Goldstone modes as func- 
tions of the strange quark mass for m u = mj = 30 MeV and 
HA 2 = 1.4 (A = 79.1 MeV): pole positions of the T-matrix at 
q = (upper panel), meson masses (center), and meson che- 
mical potentials (lower panel) . The various points indicate the 
numerical results: ty ± (circles), K + and K° (triangles), and 
K~ and K° (squares) . The lines correspond to the predictions 
from Eqs. (75), (77), and (TS). 
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For A and /„. we insert the NJL-model values obtained 
in Sects. IIV Al and IIVBI in the limit of vanishing quark 
masses. The resulting functions are indicated by the lines 
in Fig. [TUJ Obviously, they agree almost perfectly with 
the numerical calculations. Note, however, that the agree- 
ment would not be good if we had used the weak-coupling 
results for A and /„■. 

When ^ reaches the value of rrii, 



q, < 



crit \ 



mi(m 



(79) 



vanishes and meson condensation sets in. For the cho- 
sen parameters, this occurs in the kaon branch (K + and 
K°) at a critical strange quark mass mf %t ~ 145 MeV. 
For higher strange quark masses the CFL phase is no 
longer the correct ground state and the shown results 
have no physical meaning. 

Inserting the EFT expressions Eqs. ([77)) and ([75]) into 
Eq. ([79")) one finds that the critical strange quark mass 
for kaon condensation is approximately given by [26| 
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(80) 



0. For m q = 



which becomes exact in the limit m q 
this implies that kaon condensation is favored for arbi 
trarily small strange quark masses. One might expect 
that this is also the case in our model, after we found 
good agreement with the EFT predictions in Fig. [TO] 

On the other hand, this would contradict an earlier 
NJL-model study [ijj, where a nonzero critical strange 
quark mass was found, even for m q — 0. This was con- 
cluded, without explicit construction of the Goldstonc 
bosons, by studying the stability of the CFL ground state 
against partially rotating the scalar diquark condensates 
into pseudoscalar ones. Thus, with our present approach, 
we can check this result from a different perspective. 
Since in Ref. 40] indications were found that the de- 
viations from the EFT predictions are due to terms of 
higher order in the interaction, we perform our analy- 
sis using a rather strong diquark coupling, HA 2 = 1.7, 
corresponding to A = 107.9 MeV. 

In Fig. QTJ the critical strange quark mass for kaon 
condensation is displayed as a function of the third root of 
m„. The NJL-model results are indicated by the crosses. 



We also show the solution of Eq. ([79|) for the EFT masses 
and chemical potentials Eqs. ([77|) and ([78]) (solid line) 
and the approximate solution Eq. (|80[) (dashed). 

We find that the NJL points are always above the EFT 
predictions (solid line). However, while the deviations are 
small for rn q > 1 MeV, they become essential for smaller 
values of m q . In particular, we confirm that m£ does 
not vanish at m q = but goes to a finite value, which is 
about 21 MeV in our example. 
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FIG. 11: Critical strange quark mass m c " t for kaon conden- 

1/3 

sation as a function of m, , where m q = m u = m^: NJL- 
model results (crosses), EFT predictions using Eqs. (|77[1 and 
(|78p (solid line), and approximate solution Eq. (|80|l . The dot- 
ted line is based on Eq. ([5TJ with a s = 0.0203. The NJL- 
model calculations have been performed with HA 2 — 1.7, 
corresponding to A = 107.9 MeV. 



To understand this behavior, we analyze the poles of 
the T-matrix in the kaon channel for m q = as functions 
of m s . In Fig. [12] the NJL-model results (crosses) are 
compared with the EFT predictions (solid lines). As be- 
fore, the pole positions (upper panel) can be interpreted 
in terms of kaon masses (center) and chemical potentials 
(lower panel). It turns out that the latter are in fair agree- 
ment with Eq. ![78"]). On the other hand, while Eq. ([77]) 
predicts the kaon masses to be zero for m q — 0, we find 
that ma is in general non-zero and rises linearly with 
m s . As a consequence, the pole position u> K is not de- 
generate with lo k (cf. Eq. (|75[) ) and first rises with m s . 
Hence, kaon condensation does not occur at arbitrarily 
small values of m s but only for m s > 21 MeV, as already 
seen in Fig. QT] 

It should be noted, however, that, although the lin- 
ear rise of mx with m s is qualitatively different from 
Eq. ([77|) , the slope is very small on a quantitative scale. In 
our example, a linear fit mx — a, s m s yields a s = 0.0203, 
which is an order of magnitude smaller than the slope 
in Eq. {77]), y/AA/ f 2 = 0.347. Moreover, by varying the 
coupling strength, we found numerically that a s depends 
quadratically on A. This is consistent with our expec- 
tation that the effect corresponds to a higher-order cor- 
rection in the interaction and only becomes visible when 
the leading order, Eq. l[77|) . is artificially suppressed by 
choosing very small values of m q . In fact, for any realistic 
value of m u and the correction is quite irrelevant. 

Also note that our results are somewhat complemen- 
tary to Ref. [52| |. where corrections to the effective kaon 
chemical potential have been discussed. However, it was 
found there that these corrections are of the order of 
(to^/2/iA) 2 , which is completely negligible in all of our 
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FIG. 12: Kaon (K + and K°) properties as functions of the 
strange quark mass for m u = mj = and HA 2 = 1.7 (A = 
107.9 MeV): pole positions of the T-matrix at q — (upper 
panel), kaon masses (center), and kaon chemical potentials 
(lower panel). The crosses indicate the numerical results, the 
lines correspond to the predictions from Eqs. (|75[l . (|77p . and 
1)78( 1 . Note that from Eq. ()77l) niK is expected to vanish. 



examples. 

We may finally combine Eq. (|77[) with our numerical 
findings for m q — by parameterizing 



m K = J ^m q (m q + m s ) + a 2 s m 2 s 



(81) 



Equating this with the kaon chemical potential, Eq. (|78|) . 
we obtain the dotted line in Fig.[Tl]for the critical strange 
quark mass. It is obviously in good agreement with our 
numerical results. The tiny deviations for small values of 
m q can be explained by the fact that the kaon chemical 
potential in Fig. [T^] is slightly overestimated by Eq. l[75]). 
A similar effect could also play a role at larger m q , but 
there might be other terms as well. 



V. SUMMARY AND CONCLUSIONS 

We have studied the properties of pseudoscalar Gold- 
stone boson excitations in the color-flavor-locked phase 
within an NJL-type model. To that end we solved the 
Bethe-Salpeter equation in RPA. So far we have only in- 
cluded quark-quark interactions such that our Goldstonc 
boson states are in fact superpositions of diquark and 
di-hole states. 

Our results are consistent with the model independent 
predictions of the low-energy effective theory, i.e., with 
those predictions which only depend on the symmetry 
breaking pattern. We found, however, deviations from 
the values for the constants appearing in the LEET, as 
for instance the pion decay constant obtained in the 
weak-coupling limit. In fact, it was the main motivation 
of this paper to locate such deviations and to understand 
their origins. In several cases this could even be done ana- 
lytically. Although these model results are in general not 
universal, they may give important hints about to what 
extent the weak-coupling results can be trusted in the 
intermediate-density regime and where to expect major 
deviations. 

The weak-coupling limit for f n [23j is correctly repro- 
duced in the chiral limit in zeroth order in the gap pa- 
rameter A. This must be the case since this result is 
universal, i.e., independent of the specific choice of the 
interaction. For A not being small we found deviations, 
typically of the order of a few percent. We have shown 
that this is an effect of higher order in the gap param- 
eter. In this context, the momentum dependence of the 
dressed vertex function plays an interesting role. 

Next, we discussed the masses of the Goldstone bosons 
in the limit of equal quark masses. In agreement with 
the LEET prediction, we found that the meson masses 
behave linearly in the quark masses. However, the cor- 
responding coefficient A does not agree with the weak- 
coupling result obtained in HDET 23]. In the limit 
A — > 0, this should be viewed as an artifact of the NJL 
model. Probably more relevant are the deviations at large 
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values of A where we found the meson masses to be con- 
siderably smaller than predicted by the weak-coupling 
formula. 

Finally, we have studied the case of unequal quark 
masses. In general, we found a very good agreement with 
the LEET prediction for the meson masses and effective 
chemical potentials, Eqs. (|77| and (|78| . In particular, we 
found rax < m^, i.e., an inverse meson mass ordering as 
predicted in Ref. [23[ . We also confirmed that the strange 
quark mass acts as an effective strangeness chemical po- 
tential, eventually leading to kaon condensation at suffi- 
ciently large values of m s Refs. [H, H3|- However, even 
quantitatively, our model results are in an almost per- 
fect agreement with Eqs. (|77|) and (|78|) if the constants 
A and f n entering Eq. (|77p are taken from our preceding 
NJL-modcl studies in the limit of vanishing quark masses. 
Since these values do not agree with the weak-coupling 
limit, as discussed above, our meson masses are in gene- 
ral smaller than those obtained with the weak-coupling 
coefficients. 

In the limit of vanishing light quark masses, we found a 
qualitative difference. The critical strange quark mass for 
the onset of kaon condensation does not vanish with the 
cubic root of light quark masses but attains a nonzero 
value. We identified this numerically as a higher order 
effect on the kaon mass, adding a A 2 m s dependence at 
low light quark masses. 

This paper should be seen as the basis for further 
studies of the Goldstone boson dynamics in cold dense 
quark matter at non-asymptotic densities where devia- 
tions from the weak-coupling limit become visible. As al- 
ready mentioned, our simple Lagrangian does not include 
quark- antiquark interactions, which, although subdomi- 
nant, can give important corrections to the pseudoscalar 
meson masses. Moreover, we miss instanton effects, which 
have been shown to be important [H, [H| • The ultimate 
goal is to include the back reaction of the Goldstone 
bosons on the phase structure of color superconducting 
quark matter [541 ]. In the intermediate density regime, 
the ratio of the gap parameter and the Fermi energy is 
of the order of 0.25, such that the Goldstone boson exci- 
tations can have a significant effect on the ground state 
properties. 
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APPENDIX A: GAP EQUATIONS 

The Dyson equation shown in Fig. [1] reads 

S(p) = S (p) + S (p)±S(p). (Al) 



Solving for the self-energy, one obtains 

t = V(p) - s-\ P ), 

where 



A7° 





A7° 



m 



(A2) 



(A3) 



is the inverse bare quark propagator, while S 1 is the 
inverse dressed propagator defined in Eq. (|16p . Thus, 



£ = 







\A=2,5,7 



- J2 A a j 5 t a X a \ 
A ^ 5J i (Al, 



On the other hand, £ can be evaluated diagrammatically. 
In Hartree approximation, it corresponds to the quark 
loop in Fig. [T] and is given by 



S = HHT, 



d A k 1 



:2tt 



4 -Tr [TiS(k) 



(A5) 



Comparing the two expressions for S, we can read off the 
following gap equations: 



AH 



d A k 1. 



(2tt) 4 2 



-Tr 



p*l 



AA 



S(k) 



(A6) 



We also see that the contributions of scalar vertices with 
A 7^ A' and of pseudoscalar vertices should vanish to be 
consistent. Using the explicit expression for the dressed 
propagator (see Appendix [Cj , it can be shown that this 
is indeed the case. 



APPENDIX B: CHIRAL WARD-TAKAHASHI 
IDENTITY 

In this appendix, we demonstrate that the dressed ver- 
tex functions and the dressed quark propagator are con- 
sistent with the chiral Ward-Takahashi identity (%WTI) 
in the sense that, if the xWTI holds, we recover the gap 
equation. As in Sect. IIII| we restrict ourselves to the case 
of equal quark masses. 

From Eqs. ((44j) and ((47]) we obtain 



?M r 5?(i>; < ?) + 2mir 5,j{p;q) 



AiH T 



NG 

d 4 k 1, 



-Tr 



tiS(k + q) 



(2tt) 4 2 

X (q^(k;q) + 2miT 5J (k; qj) S(k) 



(Bl) 
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Hence, imposing the xWTI, Eq. ([52)1 , one gets 
S~ 1 (p + q) (75*i)jVG + (75 tj) NG S^ 1 (p). 

ffi {j5tj)NG 
+ hst^NGfi) S(k) 



Am r, 



NG 

d 4 fc 1, 
(2^)4 2" 



-Tr 



(B2) 

Using Eq. (|16[) with equal masses and gap parameters, 
one finds that the diagonal Nambu-Gorkov components 
are equal to the first term on the r.h.s., and one is left 
with 







'A* £ {TAtj+tjT A )*A 
\ ,4=2,5,7 



tjT A )X A \ 

-4=2,5,7 



(fi (75*i)iVG 

+ (75^)^^)5(^)1. (B3) 



d 4 fc 1. 



(2tt) 4 2 



-Tr 



Next we compute the sums over A on the l.h.s. for any 
of the flavor operators tj. For instance, for tj = t Tr + = 
we find that the l.h.s. is equal to 73 (A + 

A* Tyg). We conclude that the r.h.s. must vanish for all 



Tj, except for I\ = r| 7 and E; = I^g. In the first case, we 
have f , (j 5 t^+) NG + (75*7r+)wGfi = 
by comparison with the l.h.s. we obtain 



if 7V ^US 



A = 4H 



d 4 k 1 



(2tt) 4 2 



Tr 



T s 7 \S(k) 



(B4) 



in agreement with one of the gap equations (|A6[) for the 
case of equal masses. The two other equations can be de- 



rived analogously, if we evaluate Eq. (|B3[) for tj = t K + 
or tj = t^o. Moreover, the fact that most Ej must not 
contribute to the r.h.s. for a given tj can be used to show 
that scalar operators with A ^ A' and pseudoscalar op- 
erators do not contribute to the gap equation. 



APPENDIX C: DRESSED QUARK 
PROPAGATOR 

1. General case 

The dressed quark propagator S(p) is the inverse of the 
inverse quark propagator, defined in Eq. (1161) . Following 
standard methods (see, e.g., Refs. [47l |48|. l49l. [50j). we 
write 



[ (p) = S~ l (p ,p) = 7 °(p°-A(p)) ; 



(CI) 



where A(p) is a hermitian 72 x 72 matrix, which does not 
depend on p°. Thus A can always be diagonalized, i.e., 
we can find a unitary matrix U(p), so that 



A(p) = U\p)D{p)U(p), 



with 



D{p) 



V e 72 (p). 



(C2) 



(C3) 



being a diagonal matrix with eigenvalues e\,...et2- It 
can be shown that all eigenvalues are two-fold degener- 
ate, and for each eigenvalue £j, there is a counterpart — Si 
in the spectrum. This means, there are basically 18 inde- 
pendent eigenvalues. Moreover, part of the diagonaliza- 
tion is trivial because the matrix A can be brought into 
block diagonal form by reordering of lines and columns. 
The remaining blocks are in general diagonalized numer- 
ically. 

Combining Eqs. (|C1[) and (|C3[) . the propagator is fi- 
nally given by 



S(p) = U\p) 



( I 

p°-£l(p) 



\ 



U(p)j°. 



p u — £72 



mJ 



(C4) 



2. Equal quark masses 



In the limit of an exact SU (3) symmetry, we can give 
a closed expression for the quark propagator. Straight 
forward inversion of Eq. (|16|) for m u — — m a = m 
and A22 = A55 = A77 = A yields 



5' 



with 



and 



S21 



S21 S22 



75 ^ t a Xa Sn 

,4=2,5,7 



(C5) 



(C6) 



IAI 



t A t A > \ A x A ' 



m 

A , A '=2,5,7 



where we have introduced the notations 

p ± =p±Hj°, x±=p\-m 2 . 



(C7) 



(C8) 



S22 and S12 are obtained from Sn and S21, respectively, 
under the exchange /i «-> — fi and A «-> — A*. 
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The matrices Sij are 36x36 matrices representing the 
normal (i — j) and anomalous (i ^ j) Nambu-Gorkov 
components of S. Su can explicitly be written as 

1 8 

S U = S- + -(T^-S-^TaXa, (C9) 



a=0 



with 



c _ x ± (0 T + m) - |A| 2 (/ ± +m) 
* " [Pl-E^M-ED ' (C10) 

corresponding to the eigenvalue A of the gap matrix, and 

x± (j/> T + m) - 4|A| 2 (ji ± + m) 



(pg-£JT a )(p8-^ 



(Cll) 



corresponding to the eigenvalue 2 A of the gap matrix. 
The octet and singlet dispersion relations for particles 
(— ) and antiparticles (+) are given by 



= V(Vp 2 + ™ 2 tm) 2 + |a| 



(C12) 



and 



E? = V(v^+^TM) 2 +4|A| 2 , (C13) 
respectively. 



APPENDIX D: PION DECAY CONSTANT IN 
THE CHIRAL LIMIT 

In this appendix we derive a semi-analytical expression 
for the pion decay constant in the chiral limit, which is 
used in Sect. IIV Al to discuss the deviations from the weak 
coupling limit. The rf decay constant can be obtained in 
a similar way but we do not discuss this here. 

Starting point is Eq. (|43)) for the (time-like) decay 
constant in the case of equal quark masses. In the chi- 
ral limit, we have to evaluate this formula at to, = 0, 
Since both sides of this equation vanish in this limit, this 
means that the decay constant is given by the deriva- 
tive of the r.h.s. with respect to go- Moreover, we can 
employ the "Goldberger-Treiman relation" , Eq. (|55[) , to 
eliminate the coupling constant gi. We then obtain 



f: 



\dq 



(2^ 

rTr[r^)S(fc + g)(7%5iiWS(A0] 



9=0 

(Dl) 



Here we have indicated explicitly that the vertex func- 
tion depends on the momentum (see Eq. (|60p ) and is 
therefore subject to the derivative. We may thus write 



ft = ft + 5fl 



(D2) 



where 

!! = 

-A 



(2tt) 4 2 [ tW dq Q 



9=0 



(D3) 



corresponds to the contribution where the derivative acts 
on the propagator, while 



sf! = 



A / *1* I T r[^M 
(2tt) 4 2 1 dq 



9=0 



S(k) (7°76*i)wGS(*)] 



(D4) 



corresponds to the contribution where the derivative acts 
on the vertex function. 

The evaluation of ff is tedious, but straight forward. 
Inserting r' ( (0) = r- (0) from Eqs. (|57 |) or (|58 f . respec- 
tively, ti from Table [J as well as the expressions for the 
quark propagator given in appendix IC 21 in the chiral 
limit, the result for the octet mesons ("pions") reads 



f 

J 7T 



2167T 2 

2 



{^2 i°t - a+)(x + l)(3x 2 - 2x + 1) 

+ (a^ - a.g)(x - l)(3x 2 + 2x + 1) 
- a+(9x - 31) - a^(9x + 31) 



24 



ln(a| + x + 1) + ln(a 1 +x-l) 



— ln(a^ + x + 1) — hi(a 8 + x — 1) 



45 



-481n2 - 54y 2 lny 2 



My 2 
45y 2 



ln(a^~ 
x- 1 



- X 
x ■ 



-ln(c 
1 

1 



x- r 



}■ 



(D5) 



where the 3-momentum integral has been regularized by 
a cutoff A (as in the numerical calculations) and we have 
introduced the abbreviations 



af = ^(x±l) 2 +4y 2 , af 



and 



^{x±l) 2 + y 2 , (D6) 



P = Vi + y 2 - 



(D7) 



A as defined in Eq. ([55]). 

ire difficult because we i 
to know the derivative of the vertex function. For the 



Moreover, x — — and y 

The evaluation of <5/ 2 is more difficult because we need 
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flavored mesons this is encoded in the mixing angle <p(q), 
cf. Eq. ([SO)) . We can therefore write 



dq Q 



dip(q) 



g=0 



dq 



g=0 



(D8) 



where 



dip 



(D9) 



is the vertex of the orthogonal state with the same quan- 
tum numbers as the meson i (e.g., T'J^ = — 771 (rff + 

r?£)). Hence 



d i k 1 



(2tt) 4 2 



IlV[r^(fe)( 7 % 5 ^;v G S(AO]^ 



dg 



q=0 

(D10) 



APPENDIX E: EXACT FORMULA FOR A 



A 



The exact result of Eq. (|64|) is given by 
A 2 



384tt 2 



I - 481ny 2 - 321n2 + 4(o^ - a+) +44(a^ - a 

+ (16 + 3y 2 ) ln(af + x - 1) + ln(a+ + x + 1) 

+ (32 - 3y 2 ) ln(ag + a: - 1) + ln(a+ + x + 1) 
- 6y 2 In 2 

ln[4(x + 1) + a+p 2 p 3 - (3x - 13)y 2 ] 

- ln[4(a; - 1) - ar/3 2 /3 3 - (3a; + 13)y 2 ] 
+ ln[4(a; + 1) + a+/3 2 /3 3 + (3a; + 7)y 2 } 

- ln[4(a; - 1) - a 3 (3 2 (3 3 + (3a; - 7)y 2 ] ] } , 

(El) 



with a" 1 as defined in Eq. (|D6|) . 



& = V4 + Z/ 2 , ft = + V 



(E2) 



and a; = ^ and y = — , as defined in Eq. (I66|) 



Evaluating this part for the octet mesons, one finally ob- 
tains 



5fl 



1087T 2 

- 108y 2 \ny 2 - 48y 2 In 2 

+ 2(af - a£)(x 2 -x + 1) 

— 2(a^ — a$)(x 2 + x + 1) 

+ 4y 2 (a+ -a^)- lQy 2 (a+ - ci[) 



12 y l 



Aby z 



21n(a7/ + x - 1) + 2 In (a+ 
71n(ag~ + x — 1) + 7 In (a^~ 



4a; 



4a; 



90y 4 (^-— 

«o (To > 



d<p(q) 
dq 



q=0 



(Dll) 



One could also try to derive an analytical expression for 
the derivative |n=o from the polarization-loop ma- 
trix Eq. (|28p. but this is beyond the scope of this paper. 
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